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We derive a theory for transmission through disordered finite superlattices in which the interface 
roughness scattering is treated by disorder averaging. This procedure permits efficient calculation 
of the transmission through samples with large cross-sections. These calculations can be performed 
utilizing either the Keldysh or the Landauer-Buttiker transmission formalisms, both of which yield 
identical equations. For energies close to the lowest miniband, we demonstrate the accuracy of 
the computationally efficient Wannier-function approximation. Our calculations indicate that the 
transmission is strongly affected by interface roughness and that information about scale and size 
of the imperfections can be obtained from transmission data. 
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I. INTRODUCTION 



Semiconductor superlattices continue to attract sub- 
stantial interest both among fundamental and applied re- 
searchers. One motivating factor is the possibility of tai- 
loring the miniband structure [Q-^ for device purposes. 
Furthermore, a large variety of other physical phenomena 
such as the formation of Wannier- Stark ladders ^|, neg- 
ative differential conductance |^ , and Bloch oscil ations 
1^ can be observed in superlattices. 

The presence of minibands has been probed directly 
by investigating the transmission of ballistic electrons 
through short semiconductor superlattices 0,^. In re- 
cent experiments the quenching of the miniband struc- 
ture by an applied electric field was also demonstrated 
1^. Comparison of further experiments with theoretical 
calculations indicated a strong influence of scattering on 
the transmission, and it was argued that interface rough- 
ness might cause significant deviations from pure ballistic 
transmission through the sample [ p^|JTl]] . 

A good understanding of the transmission character- 
istics through short superlattices is important as these 
structures are used as energy filters. For example, in 
quantum cascade lasers, superlattice filters are used to 
selectively populate the upper energy level of the active 
region (l2|. The most straightforward way to calculate 
the transmission through a superlattice is the transfer 
matrix method |]l3|| . Alternatively, the Schrodinger equa- 
tion of the superlattice can be solved directly. These 
methods typically assume homogeneity in the direction 
perpendicular [the (x, ?/)-plane] to the superlattice. The 
momenta in the (x, ?/)-plane are then good quantum num- 



bers and decouple from the superlattice direction, reduc- 
ing the problem to a one-dimensional calculation; see, 
e.g., Ref. The one-dimensional calculation can han- 
dle fluctuations in the well or barrier thickness. How- 
ever, real samples also exhibit a lack of periodicity in the 
(x, ?/)-plane due to the presence of impurities and inter- 
face roughness. This can change the transport proper- 
ties essentially, as the states with different parallel mo- 
menta couple to each other. This (x, ?;)-plane inhomo- 
geneity can be tackled by solving the Schrodinger equa- 
tion on a mesh for the full three dimensional structure 
JlGf . Alternatively, the method of Green- functions, based 
on Ref. [0 , may be used (see Ref. for an easily acces- 
sible presentation of the method) . Recently such an ap- 
proach has been presented for a full calculation of the cur- 
rent through a resonant tunneling diode where both in- 
terface roughness and phonon scattering have been taken 
into account [Q. However, these simulations use a fine 
grid and are hence unsuitable for longer structures such 
as superlattices consisting of many wells, since the num- 
ber of grid points increases dramatically. 

In this paper we propose a new method for such cal- 
culations which significantly reduces the computational 
complexity. We treat the inhomogcneity in the {x, y)- 
plane by averaging over disorder configurations. The 
number of grid points in the z-direction is reduced signifi- 
cantly by restricting to the basis set to the Wannier func- 
tions localized in the wells. This Wannier approximation 
is shown to be valid near the resonance condition if the 
energy gap between the minibands is large compared to 
the bias and the widths of the minibands themselves. We 
compare our results to calculations on a finite grid and 
find good agreement. Our method has the advantage 
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that it corresponds to infinitely large cross-sections and 
hence, unlike the finite-grid calculations, does not show 
configuration dependent fluctuations. 

The paper is organized as follows. We present the gen- 
eral model within which our calculations are performed in 
section In section II] , we describe the approximations 
which allow us to perform practical calculations of ex- 
tended superlattice structures. Our results are presented 
in section IV and we conclude with a summary. Appendix 
shows the equivalence of the Landauer-Biittiker trans- 
mission formalism with the approach by nonequilibrium 
Green functions for the case of impurity averaging. In 
appendix ^ we justify the approximations used in sec- 
tion III. As many different symbols appear in this paper. 







for easy reference we display the frequently used ones in 
Table 1. 



II. GENERAL FORMALISM 



In this paper, we study transport through a superlat- 
tice contacted to external voltage sources via leads. We 
model the superlattice as an active central region cou- 
pled to noninteracting lead regions. This is the general 
approach described in Refs. 0,^,|o|j2| and [||. In this 
section we briefly review this approach, introduce our no- 
tation, and discuss the issue of impurity averaging. 

We divide the sample into a central region C and lead 
regions, indexed by £. The Hamiltonian is 



G-^(t,0=T*({c.W,ct(0}) 0{±{t-t')) 



(2a) 
(2b) 



Here cj (q) are fcrmion creation (annihilation) operators 
of states (j)c,i in the central region, and {• • • , • • • } denote 
anticommutators. 

In the following we consider time-independent prob- 
lems, so that G only depends on t ~ t' , and we work in 
energy-space by Fourier transforming G with respect to 
t — t' . The net current from mode a in lead £ into the 
structure is given by pc|,|2^ 



Jta — 2- 



dE 
2^ 



Tr{rto(£;) [G<{E) 

+fUE) {G'-'{E)-G-^^iE))] 



Here fga [E) gives the occupation of a state with energy 
E in lead I for the mode a, e < is the charge of the 
electron, and F/.q is a parameter describing the coupling 
between the states in the central region and the leads [see 
Eq. (^) below]. The factor of 2 is for spin. 

To describe transmission through the superlattice, we 
need to obtain expressions for the right hand side of Eq. 
(^). We do so as follows. We flrst deflne i?c,o and 
H'q, as the ordered, solvable part and the disordered 
part of the central region Hamiltonian, respectively, and 
He = Hc.o + H'fj. The retarded Green function for the 
structure is determined by the equation (see Ref. [ p2| , 
chapter 12) 



H = Hr 



H. 



ic 



(1) 



E-M 



E 



S-*(i?) G-'(i?) = l. 



(4) 



Here He and -ff^ are the terms for the central structure 
and leads, respectively, and H^c is the coupling term 
from the center to the lead £. In this paper, we ignore 
electron-electron interactions beyond Hartree so all the 
above terms are single-particle-like. 

The central structure has states with the wave func- 
tion 0cj(r), where j is the eigenstate index. We assume 
that each lead I is disorder-free so that the eigenstates 
can be separated into transverse and longitudinal parts, 
4'iaq{'r) = X£a(r)Vg(^), where z is the spatial coordinate 
in the direction towards the central structure and r is 
a two-dimensional vector perpendicular to z. The index 
a numbers the modes within a given lead. The index 
q denotes the behavior far away from the central region 
where fq{z) ~ e*'^ is assumed. 



A. Green functions and current through structure 

The current through a structure can be determined by 
the Green function of the structure in the presence of 
coupling to the leads, given by a matrix G with matrix 
elements 



The term SJ^' is the irreducible self-energy due to i?^. In 
cases where H'^ contains interparticle interactions, SJ^* 
is often very difficult to calculate; however, for static dis- 
order, simply S'J* = H^^.. The term ^Ta{E) gives the 
self-energy contributions due to the coupling of the cen- 
tral region to lead £ and mode a, 

^Ta.iji^) = 'Y{(l)cMHjc\(l^eaq){<l>laq\Hlc\<l>C',j) gTag{E) 



2tt 



k^,\H,c\<j^c,,) 9TM (5) 



where we have taken the continuum limit — > 

Li/2'K dq {Li is the length of lead £). The factor 2 
results from the two possible values ±q for a given energy 

E„ — Ef, 



Eq, 9\%{E) - ll{E 



q j^ta -r jO+) is the frec-particlc 
Green function of the lead in absence of the central re- 
gion, Eq — h'^q^ /2m*, Vq = hq/m* and E^^a is the lateral 
energy of the mode a. Here m* is the effective electron 
mass. Note that Gf^^{E) = [GJf (£;)]*, since we have a 
time-independent system. 
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The coupling parameter Tga is defined by 

_ '2'Ll{(j)CsWlc\4'e.a q{E-Ei^)) 



hv 



q(E-Ei,. 



'x{(f'ea q(E-Et^)\Hic\(t>c.]) OiE-Eea) (6) 



where q{£) — \/2rn*£/h. 8(.t) is the Heavyside function 
with e(x) = 1 for .T > and e{x) ^ for x < 0. Note 
that T£a^ij{E) ~ for E < Eia since there are no prop- 
agating states into which the central-region states can 
tunnel. 

{E) can be obtained by the Keldysh relation ||2^ ] 



where 



G<{E) = G"'\E)T,<{E)G'''^^{E) (7) 



i:<{e)^i:<{e) + Y.^Ue)- (8) 

la 

Here, is the self-energy resulting from scattering in- 
side the structure. For a fixed disorder potential, this 
term is identically zero. The term 'S,f^[E) is the self- 
energy due to the presence of the coupling to the leads. 



— * ^eaAj (E) fla {E) 



(9) 



where we have used gf^q{E) = -2ifia{E) Im js^aV-S')}- 

The occupation function fia (E) in lead £ and is given by 
the externally imposed conditions. Usually, the leads are 
assumed to be in thermal equilibrium and hence a Fermi 
distribution with chemical potential independent of 
a, is used. In contrast, the different modes can be popu- 
lated individually by injection, as discussed later, so that 
we want to keep the full function fiaiE). 



The Landauer-Biittiker approach has been used exten- 
sively to study transmission through mesoscopic struc- 
tures, and consequently many people are familiar with 
the formalism. As the Keldysh formulation is not as 
widely known, in this subsection we demonstrate the 
equivalence of the two approaches for transport through 
a system with static disorder. 

The retarded and advanced Green functions can be ex- 
pressed in terms of F via 



G^^'lE) - G'^'^^iE) ^ -iG''''{E)T{E)G''''^{E) 

(10) 

where the total scattering rate T has two contributions 



T{E) 



\E) 



-ladv 



resulting from scattering inside the structure and transi- 
tions into the leads. If the scattering within the structure 
itself is purely elastic and is treated in a particular fixed 
configuration as a potential in Eq. (|), then Eg' = 
and T,'^{E) — 0; hence we may insert Eqs. (|^) and ( p^ ) 
into Eq. yj^ and find the Landauer-Biittiker expression 

po|,pl 



■ha 



2^ / ^E^^"-^'/3(^) [hAE)~h'p{E)] 
•' e'fs (12) 



(factor of 2 for spin) with the transmission matrix 

Tea^i'piE) = Tr{rf^{E)G'-\E)r,,0{E)G-^''{E)} . 

(13) 

(There are several alternate ways to derive this result; 
e.g., Ref. [|l8| uses spatial discretization.) Note that 
Eq. ( |l3| ) does not hold if the scattering process is inelas- 
tic or the elastic scattering by static disorder is described 
by a self-energy obtained by configuration averaging. In 
both cases Sp^E) ^ in contrast to the assumption 
leading to Eq. (12). 



C. Impurity Averaging 



Eq. (|1^) is exact for a given configuration of impu- 
rities and roughness, i.e., for a specific H'^. However, 
obtaining the transmission by simulating individual con- 
figurations is not computationally efficient, and hence 
it is advantageous to average over impurity configura- 
tions. In particular, such a procedure reestablishes sym- 
metries which are broken by specific impurity configura- 
tions, thus simplifying the calculation significantly. 

After impurity averaging, we obtain 



B. Relation to the Landauer-Biittiker approach 



(11) 



G'^^(i?) ^[E + zO+ - Hc.o - ^c iE) - E ^T\E)r\ 

i (14) 

where the overlines indicate averages over disorder con- 
figurations in the central region. Note that the disor- 
der averaging introduces non-zero self-energies 'Sq (E) 
and 'Ep{E). As Sp(i?) ^ one cannot simply use 
Eq. ( p^ ) with the G's replaced by G. In order to de- 
scribe configuration-averaged elastic scattering within 
the transmission formalism, the averaging must be per- 
formed for the total transmission matrix in Eq. (13), and 
not just on the individual G"'°' and G"*^^. This proce- 
dure is similar to the calculation of bulk conductivities 
using the Kubo formula, where it is crucial to include ver- 
tex corrections which fulfill the Ward-identity (see, e.g., 
Ref. [^). We perform such a calculation in appendix ^ 
for the superlattice structure discussed in section III . We 
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use the self-consistent Born-approximation for the scat- 
tering and therefore the appropriate vertex function is 
the so-caned ladder approximation. 

The application of the more general Keldysh approach 
to calculate the current in the configuration averaged case 
is more straightforward, in that one can replace the G by 
G in Eq. (H). Therefore, in order to evaluate the cur- 

rent, we need G (E) and G (E). The general iterative 
procedure for computing these is as follows. First the 
self-energies Ti^^ and T,f (E) due to the leads are 
evaluated by Eqs. (||) and (g). As these terms are in- 
dependent of disorder configuration, these need only be 
evaluated once and then stored. With ' initially set 

equal to zero, G^ and G™ are calculated. These G's are 
used to calculate the Sc's, via an appropriate approxi- 
mation scheme. The updated Sc's are used to generate 
new G's via Eqs. and (p^, and the process is iter- 
ated until convergence is achieved. Finally, the current 
is evaluated with Eq. (|). 

In Appendix ^ we show explicitly that the ladder ap- 
proximation for the vertex function in the transmission 
formulation yields the same equations as the Keldysh 
approach within the self-consistent Born approximation, 
demonstrating the equivalence of the two methods for 
impurity scattering. Nevertheless the Keldysh approach 
seems to be conceptually easier as there is only one place 
within this formulation where an approximation is made; 
i.e., in the self-energy. In contrast, with the transmission 
formalism, errors can occur if the vertex function does 
not fulfill the Ward identity, providing a pitfall to trap 
the uninitiated and unwary. 



III. APPLICATION TO A SUPERLATTICE 
STRUCTURE 

Let us consider the superlattice structure sketched 
schematically in Fig. The superlattice consists of N 
identical wells embedded in iV -t- 1 barriers. A bias U is 
applied to the structure yielding constant potentials Ul 
and — Ul + ell at the left and right contact, respec- 
tively. In order to perform calculations we now specify 
the basis states (j)c,j{'i^) and X£a(r) for our superlattice 
structure. The lead index £ takes two different values L 
and R, for the left and right contact region, respectively. 
For superlattices with a large cross section A it is natural 
to use a basis of plane waves e*'' "'/ VA for the transverse 
coordinates (x, y) both in the lead regions and in the 
superlattice itself. Then the index a of the states in the 
leads is replaced by k and we have i?( l /_R,.k) — Ek + UL/u, 
where Ek = h^k'^ /2m. 



A. Wannier approximation for a superlattice 



tractable, we restrict ourselves to a subset of the basis 
functions of total Hilbert space, ignoring irrelevant high 
energy states. With respect to the ^-direction inside the 
superlattice we use a basis of Wannier-functions 
{n — 1 , . . . N) from the lowest miniband which are max- 
imally localized in well n ||2^. Such a basis has been 
successfully applied to superlattice transport |^^. This 
approximation, which we call the Wannier approxima- 
tion (WA), neglects higher mini-bands, and its validity 
is discussed in the Appendix There we demonstrate 
that this approximation gives good results for the trans- 
mission probability provided that the miniband width is 
smaller than the energy of the center of the miniband 
and the energy range of interest is sufficiently below the 
levels corresponding to the higher miniband states. The 
states (f)c,j within the superlattice are, within the WA, 
given by products ^'„(z)e*'^ '"/\/A, which can be labeled 
by (n,k). 

Within the superlattice, the Green-function is deter- 
mined by Eq. (El) which in the WA basis reads 



[{E -E'^-Ek- Un)Si,.w'5, 



-Tl'Jk.k' {Sn,n' + 1 + Sn,n'-l) - ^ ^Ta; nk, n'k' (-E-) 

ia 

G'n'k',mki(-^) = ^k,ki<5n,: 



(15) 



Here, Un denotes the potential in the well n (see Fig. |l|) 
which is due to an external bias. (The mean-field po- 
tential induced by the carriers in the structure can be 
added as well.) Ti is the coupling between the wells 
and E"^ is the level energy of the Wannier state rela- 
tive to the bottom of the well. For a given structure, 
we calculate Ti and i?" as follows. We consider first 
an infinite superlattice of the same composition. The 
eigenstates in the infinite superlattice are Bloch functions 
with the miniband dispersion E'^(q). E"" is then identi- 
fied as the center of the miniband d/{2Tr) J dqE°'{q) and 
Ti — d/(27r) / dqE°- (q) cos{qd) , where d is the period of 
the superlattice i.e., \Ti\ is about a fourth of the 
miniband width. Finally, Hnk n'is.' is the disorder scat- 
tering matrix element. 

If we average over disorder configurations the transla- 
tional invariance in the (x, ?/)-plane is restored, and con- 
sequently, all impurity-averaged quantities are diagonal 
in k parallel to the {x, ?/)-plane. Therefore we are able to 
use the notation GnU,mk{E) = G„„i(k, i?) and matrices 
G(k, i?) and S(k, i?) have the components G„m(k, i?) 
and I]„m(k, i?), respectively. 



Let us now consider the central region; i.e., the super- 
lattice structure itself. In order to make a calculation 
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B. Estimating the coupling and wide band limit 



C. Interface Roughness 



The coupling with the mode k in the left contact yields, 
from Eq. (|5|), the self energy 



^Lk; nki,«'k2(^) - '5n4'5n'a'5ki.k<5k2.k 

l_ r 2Ll\V,? 1 

where Vq = {ip^ {z)\Hlc\'^i{z)) is the z-dependent part 
of the matrix element for the coupling to the leads. Here 
we neglect the coupling to the inner wells (n ^ 1), 
which should be small. The right contact gives the 
same term except with replacements 6n.i Sn,N and 
Ul,Ll — ^ Uii, Lr. 

If the transmission function is strongly determined by 
resonances, only a small energy range of E k, E"" + Ul + 
Ek contributes to the transmission. In this range we ne- 
glect the q dependence of the coupling and extend the 
lower limit of the integration in Eq. ( pj| ) to — oo. Then 
we obtain for the left lead 



^Lk; nki^n'.kal^) " l5n,l<5„',l(5ki,k(5k2,k-^rL 



with 



Tr = 



hv. 



(17) 



(18) 



This approximation is often referred to as wide band 
limit. Note that this limit becomes problematic if the 
voltage drop across the first barrier becomes large, as 
this changes the relevant values of E and it cannot be 
regarded as constant [see also Appendix ^ . 

Now we want to estimate the value of |Vqp. For 
Eq Ri E°- the wavefunction fq{z) in the left lead be- 
haves like the Wannier function ^'0(2) = ^'i(z + c?) which 
is localized in a fictitious additional well on the left side 
of the structure. Now fq{z) is normalized to while 
the spatial extension of the Wannier function is given by 
Weft, which should be slightly larger than the well width, 
as the function penetrates into the barriers. Therefore we 
may set ipq{z) ~ •\/weff/iL^'o(^)- Then we can estimate 
the matrix element 



(19) 



yielding 



2woffTf 

hv. 



(20) 



For the right contact, F^^ is given by the same value. 



For ideal structures the potential i/„k,n'k' in Eq. (|T^) is 
zero due to the translational invariance within the (a;, y)- 
plane. However, interface fluctuations leading to well 
width fluctuations ^ri(r) in real samples break this trans- 
lational invariance. If interwell scattering and well-width 
correlations between different wells can be neglected, the 
averaged square of the scattering matrix element is given 
by [||j29| 



|i/„k+P,„'k|') - ^S{P)S: 



(21) 



where K is equal to the change of energy dE°-/dw 
per well width fluctuation and S{p) is the Fourier 
transformation of the well width correlation function 
(f„(r) ^„(r')) = /(r — r') which is assumed to be inde- 
pendent of the well index. The theory can be extended 
to accommodate interwell scattering and well-width cor- 
relations between different wells (which may result from 
a repetition of the microscopic interface structure over 
several superlattice periods) by the inclusion of the ap- 
propriate correlation functions (i?„ik-p,n;k^^n2k+p,n^k)- 
We use an isotropic exponential distribution /(r) = 
77^exp(— r/A) yielding 



27r 



(1 + {pXrf ' 



(22) 



where 77 denotes the standard deviation and A the in- 
plane correlation length of the well-width fluctuation. 
It is straightforward to implement more sophisticated 
distribution functions, which might be obtained from 
Monte-Carlo simulations of the growth conditions (see, 
e.g., Ref. ||3|]) or X-ray characterizations of the super- 
lattice structure (see, e.g., Ref. ||3^). Within the self- 
consistent Born approximation we obtain the self energy 



k' 



k',k 



(23) 



which provides the functional needed in the procedure 



scetched in section 11 C 



IV. RESULTS 

Let us consider the transmission of ballistic electrons 
through the superlattices considered in recent experi- 
ments by Rauch et al. The structure consists of 
N wells of 6.5 nm GaAs and -I- 1 barriers of 2.5 nm 
Alo.3Gao.7As. We obtain the band parameters = 54.5 
meV, Tl — —5.84 meV, K = 13.25 meV/nm and use 
Wcff = 10.7 nm, where we obtained the best agreement 
with "exact" calculations; see appendix |b| This value is 
somewhat larger than the well width in good agreement 
with the discussion in Section III B. We assume thickness 



fluctuations of half a monolayer 77 = 0.14 nm around the 
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nominal value and a correlation length A = 5 nm, unless 
otherwise stated. 

Motivated by the relatively sharp electron distribution 
injected into the structure, we assume that the electrons 
occupy the mode k = of the left contact at an energy 
E = Bin] i.e., we have fLkiE) = SkfiS{E - Ein) and 
fRk{E) = 0. The total current through the right contact 
is then given by 

Jfl = V J7?k = -4 E /^^ Tr {-iTr^{E) G<(E)] . 
k k^ (24) 



This can be expressed via Eq. (^^ by 



(25) 



For illustrative purpose we calculate the effective trans- 
mission T{Ein) — —JRirh/e in the following. Regard- 
ing the applied bias we assume a homogeneous electric 
field F inside the superlattice and set Ul = 0, Un — 
-{n- l/2)eFd- eFb/2 and Ur = -NeFd-eFb = eU 
where b is the barrier width. In the experiments consid- 
ered, there is no charge accumulation inside the structure 
as there is on average less than one electron inside the 
structure at a given time. If necessary such effects can be 
easily taken into account by solving the Poisson equation 
for the electron density given by 



2ttA ^ 

k 



^ ' dE Gl^i,\^,E) . 



(26) 



In Fig. y we show the effective transmission with and 
without scattering. In both cases we find a series of 
peaks, equal to the number of quantum wells, which re- 
flect the eigenstates of the superlattice structure. For 
[/ = the peak maxima reach the value 1 for the ideal 
superlattice. The broadening of these peaks results from 
the coupling to the leads and is of the order (r^^-f ri)/7V. 
In contrast the maxima are lower and the widths are 
wider for the calculation including scattering. These ef- 
fects becomes more pronounced with increasing superlat- 
tice length as the broadening due to scattering dominates 
with respect to the lead induced broadening. 

An important quantity is the integrated transmission 
for a given potential drop V 



Ti„t([/) = y" dE,^ T{E,,,-,U), 



(27) 



where the integration is extended over the whole en- 
ergy range of the band. This quantity was measured 
m Ref. Results are shown in Fig. ||. Let us 

compare the result of the calculations with (full line) 
and without roughness (dotted line) first. Without in- 
terface roughness, the function Tint(C/) is always sym- 
metric with respect to U. This can be understood 
from the symmetry property of the transmission matrix 
Ti^^t,f3{E) = Ti,p^to.{E) (see, e.g., Ref. |l8|). For an 



ideal structure, k is conserved within the structure and 
we find according to Eq. (p5|): 

Tint(C/)- / di?i„5I r(^,k)^(L,0)(i?in;f/) 
k 

= j dEin o)^(L,0)(-Bin; U) 

(28) 



— / dEin T(^L,0)^{R,0)iEin;U) 



Now T^L,o)^{R,o){Ein;U) = T(R^o)^^L^o){Ein + eU\-U) 
due to the symmetry of the structure and so we find 

ri„t(?7) = Tint(-f/). 

This argument does not hold for a superlattice with 
interface roughness as the scattering is able to transfer 
electrons from state k = 0, where they are injected to a 
finite value of k. In this case kinetic energy E^ is trans- 
ferred to the (x, ?;)-direction and the electrons leave the 
superlattice with a lower z-component of the energy Eq. 
This opens up new channels for new processes if f7 > 0; 
see also the discussion in Ref. . Therefore the function 
Tint{U) is asymmetric with respect to the bias U as can 
be clearly seen in Figure I (fuU fine). These findings are 
in excellent agreement with recent measurements pO[ |. 

In Fig. |3| we have also shown the transmission due to 
electrons traversing the superlattice without scattering 
(dashed line). This curve is obtained by neglecting the 
term S^(i;) in Eq. (|). It can be clearly seen that this 
curve is symmetric with respect to the bias and its mag- 
nitude is decreasing with increasing sample length. 

An alternative way of calculating the transmission has 
been performed in Refs. ||ll[]. There the Green functions 
were calculated for a fixed interface potential following 
Ref. [|6| . For practical reasons the size of the samples is 
relatively small. The diamonds and crosses refer to two 
different random interface potentials as shown in Fig. ^ 
which both have approximately the same statistical fea- 
tures. The data obtained for the transmission are not 
smooth for [/ > and exhibit differences between each 
other. This indicates that significantly larger areas than 
10 X 10 or 15 X 15 grid points must be used for reliable 
calculations utilizing this method, which is not practi- 
cable. In contrast the method using impurity averaging 
presented here gives a smooth behavior which, in effect, 
averages the scattered data points obtained from the pre- 
vious calculations. 

In Fig. ^(a) we have shown the integrated transmis- 
sion for different values of the correlation length for the 
roughness distributions. In the range considered we find 
that the asymmetry increases with the correlation length 
of the interface roughness. This indicates that larger is- 
lands lead to an enhancement of scattering even if the av- 
erage coverage is identical. The reason is that scattering 
events with low momentum transfer is enhanced. Such 
scattering events dominate the transport characteristics 
of the superlattice due to the energy scales involved in 
the system. Fig. ^(b) shows the increase of the asymme- 
try with the fluctuation height. The strong dependence 
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allows for an estimation of the interface quality by ana- 
lyzing the experimental transmission data. 



V. SUMMARY AND CONCLUSIONS 

We have presented a formalism to calculate the trans- 
mission of electrons through a finite superlattice in the 
presence of scattering processes. Due to impurity averag- 
ing the results are applicable to samples with large cross- 
sections. We have also shown that reasonable results can 
be obtained by restricting the calculation to a basis of 
Wannier-functions. Within this Wannier approximation 
all couplings are well defined and can be easily calculated 
from the superlattice parameters, with the only slight 
ambiguity being the effective normalization width Wcs, 
which is typically a few nanometers larger than the well 
width. 

Although we have only presented results for interface 
roughness scattering, the formalism is easily applicable to 
other elastic scattering processes, such as impurity scat- 
tering, as well. With regard to inelastic phonon scatter- 
ing, the formalism holds as well if Langreth rules [ ^J22| 
are taken into account, which provide the more compli- 
cated functionals for the retarded and lesser self-energies, 
see also Ref. jsi). Nevertheless, one encounters the prob- 
lem that the Green functions at different energies couple 
to each other. Therefore the set of equations which has 
to be solved self-consistently becomes significantly larger. 
The inclusion of electron-electron interaction within the 
mean-field model is straightforward. 

Our results show that interface roughness gives an en- 
hancement of the electron transmission for positive biases 
applied to the superlattice. The shape of the integrated 
transmissions depends strongly on the distribution of the 
well width fluctuations and allows us to study interface 
roughness in semiconductor heterostructures. This pro- 
vides a a complementary approach to the usual method 
of characterization by luminescence spectra. 
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APPENDIX A: IMPURITY AVERAGING 

In this appendix we describe the procedure of disorder 
aver agin g for the superlattice structure discussed in Sec- 
tion 



III 



The cross section A of typical superlattices is 
large enough that the transmission from L to i? is given 
by the configurational average of impurities. For the sake 



of transparency we assume that the scattering matrix ele- 
ment is diagonal in the well index and that the impurities 
are uncorrelated between different wells. The inclusion of 
both effects is straightforward and the identities derived 
below hold in a similar way if |J7(k — k')p is generalized 
to ([/„,„. (k-k')C/„,„^(k'-k)). 



1. Transmission formulation 



The formalism described here is similar to one used 
in Ref. p5| , except that that work was concerned with 
a resonant tunneling device (i.e., one well in the struc- 
ture), and the scattering was calculated only to low- 
est order. Hence only single scattering events were 
included, whereas the formalism described here takes 
multi-scattering events into account. This is important 
for superlattices, as it is unlikely for an electron to pass 
through a relatively long structure with only one colli- 
sion. 

From Eqs. (^3|) and (|l6|), the averaged transmission 
matrix can be written 

Tm^i,u^{E) - TuMiE) H(7Vk, hk,;E) (Al) 



with 



H(nk, hk^E) 



^«k, ki (^) kl (E) J^^) ki, nk(-^)' (A2) 

where the overline denotes impurity-averaging, ii — L 
or R, and we define ni{L) = 1 and ni{R) = N. This 
problem is analogous to the well-known case of impurity 
scattering in bulk material (see e.g., Ref. [^). The im- 
purity averaging introduces a self-energy lf°^^^'^^(\i,E) 
to the Green functions and, as in the bulk case, ver- 
tex corrections due to impurity potential correlations be- 
tween G"'°' and G"*^^. Within the self-consistent Born 
approximation, in addition to the impurity contribution 
to the self-energy, one must keep the ladder diagrams in 
the vertex |^ . 

We assume that every well has the same uncorrelated 
randomly distributed concentration of impurities with 
areal density nimp, and each impurity has a Fourier trans- 
formed potential U{q). Within the self-consistent Born 
approximation the retarded self-energy is diagonal within 



the well-coordinates n,m [i.e., '^cnm 
is given by 



-ret 



and 



ik,E) 



''imp 



dk' 

(2^ 



|C/(k-k')pGl(k',£;). 



(A3) 



Higher order approximations have been used in Ref. pq ] 
for the resonant tunneling diode. The impurity-averaged 



7 



Green function is obtained from Eq. (Q), and is explicitly 
given by 



vpf — 1 

G (k, E) 



Grct 




Sc"l(k, E) + Ef (k, E)6n 1 + S5?*(k, E)6nlA^) 



where n,m = I,-- - ,N are the well coordinates. The 
ladder approximation for n(nk, Ziki;i?) yields 



n(nk, £i'ki;E) 



Gn E) Ti^-[^{E) G'„^(^^) „(k, E) (5kki 

+ E "i-P / To^lt^C' - k')pG:l(k, E) 

m=l ^ ' 



(A5) 



:adv , 



xn(mk', 4ki;i?)G„„(k,i?) 



This equation can be iterated to yield n(nk, ii\ii\E). 
Note here that £iki acts merely as a parameter. As we 
assume that particles are injected from the left side into 
the superlattice with zero transverse momentum, we need 
only calculate the case £i = L and ki = 0. 



2. Keldysh formulation 

For comparison, we give below the equations which 
arise from the Keldysh formulation of this problem, 
within the same approximations described abo ve. The 
retarded Green function is determined by Eqs. (A3) and 
(A4), as in the previous subsection. In addition we have 



Sc.„(k,i?) = 



'""P / 7^ |c/(k-k')pGl(k',£;). 



(2^)2 



(A6) 



Together with Eqs. (|^) and (||) we obtain 
E<(k,i?) =*^/,k(i?)r,u(i?)<5„„,(,) 



''imp 



E 



dk' 

(2^ 



— adv , 



if/(k-k')p c:(kj^^j 



xE<(k',£;) G_(k',£;), 

The numerical expense required to find the self-c onsi stent 
solutions for S^(k, i?) and to solve for Eq. (A5) are 
the same. The final evaluation of G„„/ (k, E) — 

EmG^m(k,£^)S^(k,£:)GlJ,(k,£;) is straightforward 
and we can identify 

i n(«k, hki-E)fi,^^{E) = Gl(k, E), (A8) 



as both sides are determined by an identical set of equa- 
tions. Similarly, one can show 



- z ^ n(nk, hk,;E) = cTniK E) - G:t(k, E) (A9) 

£iki 

by using the identity G^**,*, (k, i?) — G'^J,{\s.,E) = 
d (k, E) r„ (k, E) Gtl (k, E) where r„ (k, E) = 



i[E™*(k, E) — 1]„ (k, E)] satisfies an equation very much 
like Eq. (^). The identities Eqs. (^) and (|A9|) show 



that the Landauer-Biittiker expres sion Eq. ([13) with the 
averaged transmission matrix Eq. (Al) is identical to the 
Keldysh formulation result given by Eq. (||). This ex- 
plicitly demonstrates the equivalence of the transmission 
and the Keldysh approaches within the self-consistent 
Born approximation of the scattering. 



APPENDIX B: THE ACCURACY OF THE 
WANNIER AND WIDE-BAND 
APPROXIMATION 

In order to check the accuracy of both the Wannier 
approximation and the wide-band limit, we will compare 
our results with a different approach. Calculations in real 
space have been performed in Refs. JT6|-p^ using a fine 
spatial discretization of length a. In the limit a 0, 
these calculations in principle yield exact results. Un- 
fortunately, these approaches generate huge matrices, so 
that we restrict ourselves to a one dimensional structure 
and neglect the (cc, ?/)-direction. This refers to an ideal 
superlattice, where the z and (x, y)-directions decouple. 
Numbering the discretization points with indices i the 
total Hamiltonian is then given by 



(Bl) 



with 




—J — - ^ for nearest neighbors, (B2) 



(B3) 



where a position dependent effective mass to* has been 
included following Ref. jl^. Now we assume that the 
sample is translationally invariant in the z-direction for 
discretization points z < and i > M and the coupling 
term in these regions is Ti^i±i — t < Q. Then we may 
define the region < « < M as the structure and the 
regions i < and j > M as leads within the formalism 
given above. The solutions for j < are plane waves 
sai{q]^{E)ja) and we have 



E = Ul + 2t[cos{qLa) - 1] . 



(B4) 



For Ul > E this gives an imaginary q corresponding to 
a non-propagating mode. (For a practical calculation |i| 
should be larger than \E —Ul\; otherwise Eq. ( p^ ) does 
not represent the effective mass parabola for the leads.) 
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Cutting of the leads gives a self-energy (see chapter 3.5 
of Ref. iT3): 



^00 



teyi^{i\qL{E)a\) for E>Ul 
texp(-|Ki(i;)a|) for E <Ul 



(B5) 



where kl{E) is defined by £' = C/^ + 2t (cosh Ki(i?)a — 1) 
for E < Ul- Similar relations hold for j = M with the 
mode from the right contact R. Note that this expression 
is only valid if the coupling from the lead to the central 
region is given by the same element t as used in the dis- 
cretization of the lead itself. The self energy is added to 
the potential and wc obtain the matrix equation 



[E5,, 



Hi 



G 



3 J 



(B6) 



which can be inverted to evaluate the Green function. Fi- 
nally the transmission is given by the Fisher-Lee relation 
p^, see also chapter 3.4 of Ref. Il8| 



Tr,l = 4.t^ sin\qL{E)a\ sm\qR{E)a\ \G'^UE)\ 



(B7) 



which can be inserted in Eq. ([121). 

For comparison we consider a superlattice with N=5 
wells of 6.5 nm GaAs and 6 barriers of 2.5 nm 
Alo.3Gao.7As. Then we obtain the band parameters 
E'^ = 54.5 meV and Ti = -5.84 meV. The transmis- 
sion is shown as a function of the injection energy Em 
for two different voltages in Figs. |^(a) and (b). Here we 
use WeS = 10.7 nm, where we obtained the best agree- 
ment. This value is somewhat larger than the well width 
in good agreement with the discussion in Section III B . 
We see that the transmission function contains 5 separate 
peaks which are related to the 5 states in the superlattice. 
The agreement between the approaches is quite good. If 
a bias is applied, the WA gives too high (low) transmis- 
sions for low (high) energies. The reason is the fact that 
the transmission through a barrier increases with energy, 
which is neglected in the WA. Preliminary results indi- 
cate that the agreement can be improved significantly, 
if next-nearest neighbor couplings are included both in 
Eq. (HI) and (|l|). (In this case matrix elements like 
^lL; Ik, 2kiE) and Sik;2k,2k(-£^) ^avc to be considered 
as well.) 

The results for the integrated transmission Tint{U) 
are shown in Fig. ^(c). We find, that the WA gives 
good agreement with the discrete model for the inte- 
grated transmission. The agreement becomes even bet- 
ter if a larger barrier width is used (not shown here). In 
Fig. |^(d) we examine the length dependence of the inte- 
grated transmission calculated within the discrete model. 
We find that the function Tint{U) becomes almost inde- 
pendent of N for large N. The results from the WA 
are almost indistinguishable and not shown here. Note, 
that Tint{U) is a symmetric function with respect to the 
bias, which can be shown analytically using the symme- 
try properties of the transmission matrix, see Section IV . 



In order to estimate the range of validity of the WA, 
we have also considered different superlattices. By de- 
creasing the barrier width to 1 nm, we have generated 
a strong coupling between the wells. Here the mini- 
band width is slightly larger than the center of the mini- 
band. Best agreement between the approaches is found 
for WcS = 13.5 nm. This value is larger than the one ob- 
tained above as the Wannier functions are less localized 
due to the small barrier width. As shown in Figs, ^(e) 
and (f) the results from the WA deviate clearly from the 
"exact" result in this case. 

Finally we considered the case of a larger well width, 15 
nm. The calculated miniband width (4Ti) is about 75% 
of the the center of the miniband Nevertheless, the 
agreement between both approaches is still satisfactory 
for low biases, as shown in Fig. ||(g). The second mini- 
band extends from 44.8 meV to 84.7 meV in this case. 
Its influence can be seen in the integrated transmission. 
Fig. |(h). For \\J\ > 0.024 V the apphed bias is larger 
than the gap between the lowest and second miniband. 
Then the coupling between the bands becomes important 
and the integrated transmission increases with bias for 
\U\ > 0.03 V for the calculation in the discrete basis. 
Naturally this effect is not accounted for in the Wannier 
approximation due to the restriction to the lowest mini- 
band. 

In conclusion we find that the Wannier approximation 



toget her with the wide band limit from Sees. [II A and 
III B gives good results if the miniband width is smaller 
than the energy of the center of the miniband and the 
applied bias is smaller than the gap between the mini- 
bands. 
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(.,(' 


Lead index 


a, (3 


Subband index (in leads) 


C 


Central region 


\(t>taq) 


State in lead £, subband a with 


wavenumber q 




State in central region 




Index of central region states 


n, m 


Index of quantum well 


fea{E) 


Distribution function of lead £, subband a 


G 


Green function in central region (Matrix 




with respect to i, j or n, m) 


S 


Self-energy in central region 


r 


i(S'-et _ s*^^) 


rct,adv,< 


Retarded, advanced, "<" component of 




Keldysh functions 


G 


Impurity average (in central region) of G 


k 


Wavevector in (x, y)-plane 




[i.e., plane ] to superlattice interfaces] 


Eq 


"Longitudinal" kinetic energy in lead 




h^q^/2m (assumed i,a independent) 




"Transverse" kinetic energy of lead £, 




subband a 


E°- 


Energy of middle of miniband in central 




region 


Vg 


Longitudinal velocity of particle hq/m 




in lead 




Applied potential in lead i 




Length of lead £ 


Jia 


Current into central region from £, a 


e 


charge of the electron (e < 0) 



Table 1. Table of various symbols used in this paper. 
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FIG. 1. Sketch of the structure considered. The energy lev- 
els are indicated for Ek = 0, which has to be added for finite 
parallel momentum k. 
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FIG. 2. Transmission through superlattices for different 
lengths and biases: The full line depicts the transmission cal- 
culated from Eq. for a superlattice structure with inter- 
face roughness. The dotted line gives the transmission for an 
ideal superlattice without scattering. 
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FIG. 3. Integrated transmission from Eq. through su- 
perlattices for different lengths: The full line depicts the 
transmission calculated from Eq. for a superlattice struc- 
ture with interface roughness. The dashed line denotes the 
part of transmission without scattering. The dotted line gives 
the result for an ideal superlattice without scattering for com- 
parison. The crosses and diamonds give the results calculated 
for the same structure within the model of Ref. ^| for two 
different realizations of the interface roughness. 
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FIG. 4. Two different realizations of interface struc- 
ture used witflin tfie model of Ref. Tlie left one 
is used for the crosses and the right one for the dia- 
monds in Fig. ^. Here we show the distribution of the 
second well, the other distributions have identical statis- 
tical properties. The spatial discretization is assumed to 
be 5 nm. Both distribution give a spatial correlation 
((eW - (C(r)))(e(r') - (C(r')»> ^ 0.52exp(-|r - r'|/5nm) 
within the next 2 neighbors, (^(r) denotes the local fluctua- 
tion.) 
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FIG. 5. Integrated transmission through superlattices with 
different correlation lengths (a) and different heights (b) of the 
interface roughness. 
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FIG. 6. Transmission through ideal superlattices without 
interface roughness, (a-c) and (e-h): Comparison between 
the "exact" model with flne discretization (full line) and the 
Wannier approximation (dashed line) for superlattices with 
5 wells, (d) Comparison of the integrated transmission for 
different superlattice lengths using fine discretization. The 
superlattice has GaAs wells with widths of 6.5 nm (a-f) and 
15 nm (g-h), and Alo.3Gao.7As barriers with widths of 2.5 nm 
(a-d) and 1 nm (e-h). The calculated parameters for the cen- 
ter of the miniband i?" and the coupling T\ (about a fourth 
of the miniband width) are given on top of each graph. 
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